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Abstract 

A semilinear heat equation ut = Au + f{u) with nonnegative initial 
data in a subset of L^{Q) is considered under the assumption that / is 
nonnegative and nondecreasing and Q C R". A simple technique for 
proving existence and regularity based on the existence of supersolutions 
is presented, then a method of construction of local and global supersolu- 
tions is proposed. This approach is applied to the model case f{s) = s*", 
<f) G L''{Q): new sufficient conditions for the existence of local and global 
classical solutions are derived in the critical and subcritical range of pa- 
rameters. Some possible generalisations of the method to a broader class 
of equations are discussed. 



1 Introduction 

In this paper we consider the question of the local existence of solutions to the 
semilinear heat equation 

Ut — Au = f{u) in r^T, 

u^O on dnx{0,T), (1) 
m(0) = </> in 

where is a smooth domain in R" and fir = {0,T). We assume throughout 
that the source term / : [0, oo) i— > [0, oo) is continuous and nondecreasing. The 
initial data is always assumed to be an element of L\{n) i.e. the set of a.e. 
nonnegative, integrable functions on 57, or some subset thereof to be specified 
in the course of the presentation. 

We are interested in classical solutions, that is functions u : flT i— K with 
u(0) = (j) that satisfy the equation pointwise with derivatives understood in 
their classical sense and continuous up to the boundary dfl. It is known that 
when 4> G L°°{n) then the problem admits a local classical solution, see e.g. [3]. 
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For unbounded data however statements of such generahty are not available 
and more specific model problems are considered with initial data belonging to 
some well-undestood Banach or Hilbert space and some particular choice of the 
source term allowing for analysis. 

As indicated in the opening paragraph we also impose restrictions on the 
initial data and the source term. However apart from nonnegativity and inte- 
grability we do not impose any additional conditions on the data allowing the 
existence argument to identify the subset of admissible data. 

We follow the standard practice of tackling the problem ((T)) indirectly via 
the associated integral formulation, the so-called variation of constants formula 

u{t) = S{t)(j)+ f S{t- s)f{u{s))ds. (2) 

In this formulation {S{t)}t>o represents the heat semigroup associated with the 
domain 51 and is defined as 

(S(m{x)= I K{t,x,y)c^{y)dv for x 6 H, (3) 

where > is the heat kernel for the domain 57, see e.g. In other words 
the function v(t) — S{t)(j) solves the linear heat equation 

vt- Av ^0 in n, t> 0, 
u = on dn, t>0, 
v{0) = (p in n. 

Every solution of the integral formulation ^ that is bounded for t > must 
solve ^ . Hence the route to classical solutions of the original partial differential 
equation proceeds through analysis of bounded solutions of the variation of 
constants formula. 

The existence argument involves the operator 

-nvW) = S{t)cj, + f S{t - s)f {v{s)) ds (4) 
Jo 

which we will use to define a sequence of functions converging to a solution 
u = ^[u] on 57t. The method used relies on monotonicity of the operator 
^ which results from positivity of the heat semigroup and monotonicity of 
the source term /. It should be mentioned that a similar technique was used 
by Weissler in his investigation of the global solutions to the model problem 
/(s) = sP in the whole space, see |7j. 

The statement of existence of local classical solutions is divided into two 
parts. The abstract existence proof assumes only that we have a supersolution 
i.e. a function w such that ^[w] < w. Existence of such an object combined 
with monotonicity of ^ will turn out to give us existence of a solution of the 
integral formulation. Hence the problem of proving existence of a solution re- 
duces to finding a supersolution which, as we will show in the sequel, is relatively 
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straightforward in many interesting cases. An immediate benefit of the tech- 
nique is that the regularity of the solution is then inferred from the properties 
of the supersolution e.g. if the supersolution is bounded for positive times then 
it follows from the standard regularity results that the solution is necessarily 
classical, see [1]. Hence the regularity is "for free" and does not have to be 
obtained by means of bootstrap arguments as e.g. in [21 [J. 

In the second part of the argument we provide a practical way of finding 
suitable supersolutions. We will arrive at certain structural conditions specifying 
the initial data that allow in certain examples to define a supersolution and thus 
deduce existence of a solution of ©. In particular we will see that the central 
object is the integral 



which involves all the parameters of the problem. 

Next we test our findings on the model problem /(s) — with initial data 
in (r2) . The problem of existence of local and global classical solutions is well- 
understood and explicit criteria discriminating well and ill-posed problems are 
known and presented in literature in the form of inequalities involving critical 
exponents. For this model problem the critical exponent is given by qc = "^^^'^'^ 
and well-posedness is usually analysed in three regimes: supercritical q > qc, 
critical g = gc > 1 and subcritical q < q^. These relations reflect the balance 
between the assumed regularity of the initial data, the growth of the source 
term and the dimension of the domain. 

The existence and uniqueness results are obtained using Banach's fixed point 
argument in a suitable metric space of curves in L'^(ri) satisfying certain norm 
decay properties, see [US]. The original proof relying on the Banach contraction 
theorem appeared in [5]. In the same paper the author offered also a second 
method of proof based on the positivity properties of the heat semigroup that is 
similar to our approach. In the discussion that follows the author compares the 
two existence theorems to conclude that the positivity based argument recovers 
the existence result in the supercritical range of parameters, but does not repro- 
duce the results in the critical and subcritical cases. In particular the positivity 
method as described in [5] does not allow one to deduce the existence of local 
classical solutions for all initial data in L'^iQ) in the critical case. The proof of 
this result relies on a delicate application of smoothing estimates in Lebesgue 
spaces followed by a bootstrap argument, see [5] for details. 

The particularly appealing property of the positivity argument is that it 
yields pointwise bounds on the solution in the form of space-time profiles whereas 
the usual approach offers a rate of decay of Lebesgue norms instead. Moreover, 
as we will show later on, the positivity based existence theorem is stated in 
terms of integrability properties of the map 1 1— ^ ||S'(i)(/)||ioc(f2) rather than crit- 
ical exponents associated with smoothing estimates in Lebesgue spaces. An 
assumption of this kind is much more general as it does not presuppose the 
space of initial data. For example the existence claim in the supercritical case 
mentioned above requires only that the integral JJ" ||S'(s)0||^"^\j^Nds is bounded 




(5) 
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for some t > which may be reahsed irrespective of 4> being in any particu- 
lar Lebesgue space, whereas the critical exponent condition n{p — l)/2 < g is 
specifically related to the context of Lebesgue spaces. It is therefore desirable to 
extend the method based on positivity to the critical and subcritical range. This 
is done in Section 0] where in fact the three ranges of parameters are addressed 
in a unified approach. 

The method also provides sufficient criteria for global existence without any 
additional effort. As mentioned before, global solutions for the model case 
/(s) = where investigated in [7]. The global existence result found there relies 
on the smallness of ||S'(s)(/)||^^^^Q^ds in the supercritical case or smallness of 
II0I|l<!(si) in the critical case. Again in the supercritical case the argument uses 
the positivity of the heat semigroup and provides a space-time profile. However, 
in the critical case a norm based technique is used instead. Our construction of 
the supersolutions provides new sufficient criteria for existence of global classical 
solutions in the critical and subcritical range of parameters. 

The general construction of supersolutions contained in Proposition 13.21 is 
by no means restricted to the case of simple polynomial nonlinearity and may 
be applied to other source terms of interest. In the last section we discuss some 
possible extensions. 

2 Existence result involving supersolutions 

Let denote the set of nonnegative, almost everywhere finite, measurable 
functions on fJ^. For our purposes it suffices to define a solution to be any 
u G M.'^ such that J^[u\ = u a.e. in ^j-- Any w £ satisfying ^[w] < w 
(resp. ^[w] > w) will be called a supersolution (respectively subsolution) . 

The following result serves to reduce the problem of showing existence of a 
solution to that of showing the existence of a supersolution. 

Theorem 2.1. Assume that f : [0, oo) i— > [0,oo) is continuous, nondecreasing 
and let (j> g L^{fl). Then the operator ^ admits a solution in fix if and only 
if it admits a supersolution in fix- 

Proof. By definition every solution is at the same time a supersolution and 
we only need to show that existence of a supersolution implies existence of a 
solution. 

Our existence argument relies on positivity and monotonicity properties of 
the operator Take u,v G M^, u < v a.e., u(0) < t;(0) a.e. in £7 such that 
^[u],^[v] e Xj. Then 

^H(i) - ^Mit) = S{t)iviO) - u{0)) + f S{t - s) if (vis)) - f{uis))) > 

follows since by the monotonicity assumption on / we have f (v{s)) — f {u{s)) > 
and the heat semigroup is positivity preserving, see [1]. 
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Starting with the supersolution w we can create a sequence {wk = '^''[w]}k>o, 
wq = w. Now due to monotonicity and positivity of ^ we find that 

.^[w] <w^.^^[w]< .^[w] ^ ■■■ ^ .^''+^[w] < ^'^[w] . . . a.e. in VLt. 

Thus over almost every point (t, x) S f^T we have a nonincreasing, nonnega- 
tive sequence {'Wk{t,x)}k>o which allows us to define the candidate solution by 
taking the pointwise limit 

^00(^,2;)= lim Wk{t,x) whenever it exists. 

It remains to show that ^[woo] = Woa a-e. in fix- With the integral rep- 
resentation ^ ill mind we can apply the monotone convergence theorem to 
get 

lim ^[wk]{t,x) ^ S{t)(f) + / / K{t~s,x,y) lim f{wk{s,y))dy ds 
a.e. in 17^. By continuity 

lim /(wfe(s,y)) = / I lim Wk{s,y)] = f{woo{s,y)) a.e. in f^T- 



Thus we have continuity of ^ at Woo i.e. 



lim ,^[wk\ = ^ 



lim Wk 

/c— >-oo 



^[woo] a.e. in Qt- 



Now, due to the construction of the sequence we have Wk+i — '^[wk] and so we 
finally arrive at 

Woa = lim Wk+i = lim ^[Wk] = -^[Woa] 
a.e. in fix, which means that Woo is a solution. □ 



3 Identifying supersolutions 

Theorem 12.11 presupposes existence of a supersolution without indicating how 
to construct one. In this section we will present a practical way of finding a 
candidate supersolution. The algorithm will be then tested on a specific well 
researched model case of a polynomial nonlinearity and compared with known 
results. 

The guiding principle is that equation Q is a perturbation formula in the 
sense that the solution is assumed to be given by the linear heat flow S{t)(j) 
which is then supplemented by a correction term that accounts for the effect of 
the nonlinearity. Hence we would expect that the contribution of the integral 
term is initially small in some sense as compared to S(t)(j). The aim of this 
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section is to devise a method of modifying the integral part of the formula ^ 
so that the resulting object is a prospective supersolution. 
Suppose that there exists a solution u : fir >— >■ [0, oo), then 

u{t) ^ S{t)(f)+ f S{t-s)f{u{s))ds>S{t)(l) 
Jo 

and if we now apply the operator ^ to both sides repeatedly then due to 
monotonicity we will find that 

^''[S{-)(f\{t) < u{t) for every fc > 

i.e. every ,^'^[S{-)(j>]{t) is a subsolution. 

Our strategy now is to take the first nontrivial subsolution, which is ,^[S {■)(/)] (t), 
and by bounding it from above derive a candidate supersolution. Hence broadly 
speaking we will replace 

^[S{-)cl>]{t)^S{t)^+ f S{t~s)f{S{s)cl>)ds 
Jo 

with ^ 

w{t) = S{t)<i) + f S{t-s)F{s)ds 
Jo 

for some F e such that 



S{t~s)f{S{s)(j))ds< / S{t-s)F{s)ds 
Jo 

hoping that w defines a supersolution on VLt for some T > 0. The optimal way of 
finding a supersolution depends on / and no one general recipe seems available. 
Nevertheless there is a simple generic approach which gives satisfactory results. 

Proposition 3.1. Suppose there exists an integrable function h : [0, T] i— > [0, oo) 
and some ip G L\^(Vl) such that for t G [0, T] 

f (^S{t)(t> + S{t)iij h{s)d^ < h{t)S{t)il;, (6) 

then ^ 

w{t) = S{t)(f> + S{t)tl; [ h{s)ds 



Jo 

is a supersolution on fix- 
Proof. For t G [0,T] we have 

^[w]{t) ^ S{t)(t> + S{t- s)f (^S{s)(l) + S{s)iP I h{r)dr ] ds 
<S{t)(j)+ f S{t- s)h{s)S{s)i;ds. 
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Since h does not depend on the space variables and due to the definition of the 
heat semigroup we can write 



S{t- s)h{s)S{s)i;ds 



h{s)S{t - s)S{s)'il^ds 



h{s)S{t)2pds = S{t)tl; / h{s)d 



and so ^ 

.^[w]it) < S{t)(j) + Sit)^p [ h{s)ds 

Jo 

for t G [0, T] i.e. ^[w] < w on Ht as required. 

The next step is to propose a way of defining function h. 



□ 



Proposition 3.2. Let (p G L^{fl) and suppose that there exist constants A > 1, 
T > and functions ip G L^{il), g : [0, oo) i~> [0,oo) and h : [0,T] ^ [0,oo) 
such that 



S{t)(l) <g{S{t)i{;) am 
for t lE [0,T], where 

h{t) = 



1 



gism) 



L°°(n) 



/ h{s)ds < A 
Jo 



(7) 



fiAgiSit)^)) 



L°=(n) 



Then w(t) = S{t)(/) + S{t)ijj Jg h{s)ds is a supersolution for t G [0,r]. 

Proof. The conditions hsted above serve to ensure that the inequahty ^ holds. 
We have 



h{s)ds 



□ 



/ ( S{t)<p + S{t)^j I h{s)ds ]<f[g iS{t)^p) + S{t)^ h{s)ds 



<f{Ag{S{tm<h{t)S{t)i;. 



L°°(n) "'0 



The function g used above is intended to represent an operation that allows 
us to extract the presupposed regularity of the initial condition. In the next 
section we will see how one may use the Jensen inequality to take advantage of 
the L''- integrability of (p. More precisely Lemma informs us that for r > 1 
and (p nonnegative we have {S{t)(pY < S{t)(j)^'. Thus for cp G i'^(i^) we can 
write 

s{t)^^i{sm')^ < (Sim')' 



7 



so that the correspondence to the above is given by ■0 = and g{s) — s^/'. 

Observe that in the results above the apphcabihty of a given supersolution 
is vahd as long as the inequalities ^ and (O hold. Hence if we consider the set 
of initial conditions for which T = oo then our results turn into the statements 
of global existence of classical solutions for small data. In this case the global 
supersolution ■w{t) — S{t)4> + S{t)ilj h{s)ds may be for example used to obtain 
a pointwise bound on the asymptotic profile of the solution. Note however that 
the smalhiess of data is understood in the sense of said inequalities and may differ 
from the statements found in the literature where the smallness is sometimes 
defined in terms of the norm of the initial condition in some relevant functional 
setting, cf. [1[7]. 

4 Example - f{u) = 

In this section we restrict our attention to the case (f) € L^_(r2), q > I. For p > 1 
consider the following model problem: 

ut = Au + vP in Vl, t > 0, 

u^O on dn, (8) 
u(0) = (p in n. 

Existence and uniqueness results for this problem were established by Weissler 
and then augmented by Brezis and Cazenave |2] . For comparison purposes we 
present an abridged version of their findings concentrating only on the existence 
of classical solutions for nonnegative data. 

Theorem 4.1. Let (p G L'^(fi), q > 1, and suppose that one of the following 
cases holds: 

1. supercritical: q > n{p — l)/2, 

2. critical: q = nija — l)/2 > 1, 

3. suhcritical: n{p — l)/2p < q < n{p — l)/2 and 

1 Tl 

limp"S'(t)0||^p,(f^) =0, '^^ere a = ^— Y - — , (9) 

then there exists T — T[(j)) > such that the problem (0) has a local classical 
solution on ^It ■ 

As mentioned in the introduction the method of proof relies on the contrac- 
tion argument in a carefully chosen space of curves followed by the bootstrap 
argument for regularity. In [5] the supercritical case was also resolved using 
positivity of the heat semigroup but the critical range was not covered. Here 
we present a simple way of extending the positivity-based method to include 
critical and subcritical ranges. 
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Before we present the full result let us see how the local existence and regu- 
larity in the supercritical range n{p—l)/2 < q can be recovered in a short direct 
computation. To do this it sufRces to use the standard L'^ — V smoothing esti- 
mate [3]: 

l|5(t)^||L.(o) <i"^^''^)||0||L.(n), (10) 

where 1 < q < r < oo. 

Let A > 1 be a constant. We will show that in this case AS{t)(l) is a 
supersolution for ^ on some sufficiently small time interval. We have 

nAS{-)m < Sit)<p + S{t - s) (^A'^\\Sis)<t>rL:.\nyS{s)<p) ds 

- S{t)cb (^1 + AP \\S{s)cf>r-J^,,^ds^ < AS{t)^, 

whenever ||S'(s)0||^^^^j^jds is finite and t is small enough. In the considered 
supercritical range the smoothing estimate yields 

ii5(.)0iir4^)d. < mr.in/'"^ 

and for every (f) G there exists a time T — T (||'/'||l9(o)) such that 

i+A^ur-i,,/-^ <A (11) 

for t e [0,T]. 

It is worth mentioning that in order to obtain the global supersolution for 
small data the standard smoothing estimate is not enough. For every A > 1 and 
1 1 1 1^9(0) there will be a time such that (|lll) ceases to hold. However the global 
validity of the supersolution depends on the smallness of the \\S{s)if)\\Pj~}^^-^s^ds 
rather than the bound obtained with the smoothing estimate. The explanation 
of the discrepancy lies in the fact that the smoothing estimate works well for 
small times but is too crude for large times. When t is large one should use 
the bound involving exponential decay, see [4J p. 441. We will now leave the 
discussion of global supersolutions and come back to it after the more general 
result is presented. 

In order to include critical and subcritical cases a slightly subtler choice of 
the supersolutions is required. We need two additional results. The first one 
is a particular case of the standard smoothing estimates, see |2J for the full 
statement and the proof. 

Lemma 4.2. Given 4> G L'^{^), 1 < 9 < oo, then for every K > there exits a 
time T ~ T{(j), K) > such that 

t^\\S{t)^L^(^n)<K (12) 

forte[0,T]. 
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The second result involves the interplay between convex functions and the 
heat semigroup as found in [5] [8] . 

Lemma 4.3. Let <j) > Q he a measurable function on and r > \, then 

{S{tW < S{tW. (13) 

The main result of this section follows from a combination of the Proposition 
13.21 and the above mentioned lemmas. 

Theorem 4.4. // 

WSmn'^^^n) jl \\S{sWiSin) < (14) 

for i G [0, T], < T < oo, where Cp = , then problem ^ has a classical 

solution on fix ■ 



Proof. We simply identify the elements of the Proposition 13.21 in the current 
setting. Lemma 14.31 provides us with the explicit choice of 17(5) = s^/'' and 
%j} = (f)'i so that the supersolution is given by 



w{t) = S{t)<p + APS{t)<P'^ l|5(s)0^ II^L(a) ds 
and the condition ([T]) reads 

II^W0'IU(n)/*II^WIlA^)d^<^ 

for some A > 1. Now, the right hand side attains its maximum on (1, 00) equal 
to yielding the desired result. □ 



The supercritical case is immediately recovered once we set g = 1 in (|T] 
Let us now turn to the critical case n{p — l)/2 = q > 1. Due to the smoothing 
estimate we have 



and 

ii„-„.^<: (15, 

since q > 1. Thus we can continue the estimate to get 

/^^t-iil-?-!) r < P^zIk^ < Cp, (16) 

Jo 9—1 

where the last inequality follows from the Lemma 14.21 for times small enough. 
Hence Theorem 14.41 immediately implies Theorem 14.11 with the additional ad- 
vantage of the space-time profile given by the supersolution. 
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In the subcritical range n{p — l)/2 > g it is known that the problem ([5]) is 
ill-posed in L''{U) in the sense that for any <j) e L'(ri) one can find a sequence 
(pn of nonnegative smooth functions convergent to (j) in but such that 

blow-up times of the associated solutions converge to zero, see [T]. Thus in this 
range not every initial condition gives rise to a solution. It is therefore expected 
that the calculation (ITS)) cannot be extended beyond the critical combination 
of parameters. Nevertheless the condition does not preclude initial data 
in the subcritical range. It is clear however that the decay of ||S'(i)(/)||ioo(Q) 
as described by the smoothing estimate for a general element of L''(0) is not 
enough to ensure that the inequality in question is satisfied. Thus we are led to 
consider a subset of initial data characterised by a faster decay of the supremum 
norm under the action of the heat semigroup. The following result will enable 
us to relate our findings to the third option of the Theorem 14. II 



Proposition 4.5. Let q > 1 and suppose that (p G ^+(^) such that there 
exists T = t((/)) such that 

t^^\\s{tw\\l^^^^<cf'\ 

then problem (0) admits a local solution. 
Proof. We have 



g- 1 / p— g 



Sv- 



1 ds = C„ 



as required by the condition □ 

Observe that for a general element (p 6 L'^{ft) the Jensen inequality (|13p 
yields 

\\Sm\L..(n) < ||5(i)<^^|li(a) 
and at the same time the smoothing estimate (|10p implies that both norms 

share the same decay rate t . If we then assume that the initial condition 

(j) £ is such that 

i^ll^W^'llifo) ^0 as t^O (17) 
in place of assumption ([S]), then we have 

and so 

i^||5(0<^'^||^) <2" ||5(V2)||i(,,)^0 as t->0. 
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Thus Theorem 14 . 1 1 and Proposition |33] agree for initial data satisfying condition 
(|17[) . It would be interesting to determine whether there are initial conditions 
that comply with the requirements of Theorem 14.11 but not of Proposition 14.51 
Interestingly the sufficient condition given in the Proposition 14.51 bears a 
striking resemblance to the necessary condition that any nonnegative integral 
solution of ([S]) must satisfy, see [8j. In the article quoted we find the following 
result. 

Proposition 4.6. Suppose that u G M.^ is an integral solution of the equation 
0i in the sense of the variation of constants formula. Then necessarily 

t^\\S{tm\L^(^n) < 

fort e (0,r] and x eVl. 

This intriguing formal similarity between the sufficient and the necessary 
conditions requires more research. 

Let us conclude this section by restating the existence result in a unified 
form. The inequality (fT4|) gives a well-posedness condition of a more general 
form that the one involving critical exponents. To identify the admissible data 
we can define the operator 

p(</.,T)= sup (\\s{tm^ f wsism'^^^^^ds) 

te[o,T] \ ^ ' Jo J 

so that we can say that problem ([5]) has classical solution for every initial con- 
dition in the set 

X = e L\(Sl) : P{(\),T) < Cp for some T e (0,cx)]} , 

where the global supersolutions are guaranteed for all for which P{(j), oo) < Cp. 

With this notation we can rephrase our findings and their correspondence to 
the existing results. In the supercritical and critical range we find that L'^(i}) = 
X. Indeed, the calculation above states that for this choice of parameters every 
g Ll{n) satisfies P{(t),T) < Cp for sufficiently smaU T so that i^(l^) C X. 
On the other hand X C L'^(il) by definition. In the case n{p—l)/2q > 1 instead 
of equality we have a proper inclusion X C which reflects the fact that 

not every initial condition in i^(f2) gives rise to a solution. 

5 Final comments 

The method of supersolutions shares some features with Banach's fixed point 
argument which is a standard technique in the field of well-posedness of PDEs. 
There are however essential differences which may render the method described 
here easier to use. Observe that there is no need for constructing a particular 
space for the method to work. The fixed point argument requires a contraction 
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space, i.e. a set of curves that when equipped with some appropriate metric 
becomes a complete metric space and such that the operator ^ restricted to this 
space is a strict contraction. Finding an appropriate choice of the contraction 
space may prove difficult in many cases and may require great intuition and skill 
in manipulating inequalities. This point is illustrated by the construction of the 
contraction space for the model problem ([5]) in [S] ■ The difficulties arise partly 
because the contraction argument yields not only existence but also uniqueness 
and so the contraction space needs to accommodate properties necessary for 
proving both. The method described in this paper does not address the issue 
of uniqueness, which we leave for a dedicated discussion elsewhere, but at the 
same time asserts the existence of solutions under arguably milder conditions 
than those usually required in fixed point arguments. 

One of the most important features of the technique is that the initial data 
(\) is, to begin with, only required to be nonnegative and integrable and it is 
the procedure itself that imposes further restrictions in the course of estimates. 
Hence, it is the source term / that defines the space of admissible initial data 
rather than the data restricting the choice of the source term. In consequence 
we may find that the given problem possesses local solutions for initial data that 
forms a set that is not necessarily a linear space. Such a situation should in 
general be expected when the equation at hand is a nonlinear one, nevertheless 
at the time of constructing the contraction space one is usually interested in 
well-posedness in some "well-behaved" space like a Hilbert or a Banach space. 

The method yields pointwise bounds in the form of space-time profiles which 
may prove useful in investigations of such topics as decay rate, blow-up rate, 
blow-up profiles, blow-up sets and asymptotic behaviour in instances when 
global supersolutions are available. 

The method affords many possible generalisations. All the results hold if we 
relax the requirement that the domain is smooth and demand instead that 
it satisfies an exterior cone condition, see [4j p. 440 and references therein. In 
the analysis of the integral equation ^ the heat semigroup may be replaced 
by a more general one as long as the requirements of the theorems are satis- 
fied, namely that the family {»S'(t)}j>Q is a positivity preserving semigroup on 
L^i^i) such that S{t)(\) G for any G L^(yL). The restriction that / 

should be nondecreasing may be weakened since the necessary ingredient of the 
existence proof is the order preserving property of the operator ^ restricted 
to the sequence {u;fc}fc>o which is guaranteed if / is nondecreasing but may be 
realised without this assumption. As far as the existence proof is concerned 
the requirement that / be continuous may be weakened since a nondecreasing 
function is almost everywhere continuous and the result is meant to hold almost 
everywhere. 

In the example of the polynomial nonlinearity ([8]) the key ingredient that 
allowed us to extend the monotonicity argument to cover the critical and sub- 
critical cases was the application of the Jensen inequality . Thanks to this 
inequality it was possible to make better use of the assumed integrability of 
the initial data. Note that the general construction of supersolutions suggested 
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in Proposition 13.21 does not specify how the extra regularity is obtained. The 
Jensen inequahty is particularly apt in the context of Lebesgue spaces. It would 
be interesting to see how to extract the required information in other functional 
contexts e.g. Sobolev spaces or Orlicz spaces. 

In the subcritical range the results obtained with the supersolution method 
resemble those described in |S1[H]. At this stage however the similarity is formal 
and more research is needed in order properly to relate the two approaches. 

Lastly, it would be desirable to relax the nonnegativity assumption on the 
initial data and to guarantee uniqueness. These issues require a separate treat- 
ment and will be discussed elsewhere. 
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